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Abstract 



3; 

■ Within the heavy baryon chiral perturbation theory we derive the leading chiral cor- 

00 ! rection to the nucleon GPD at £ = 0. We discuss the difficulties of consideration non-local 

light-cone operators within heavy baryon approach and the methods to solve the difficul- 
^ cj ties. The consideration of the chiral corrections directly for non-local operators allows to 

£^■1 resolve the ambiguity of the inverse Mellin transformation. In particular, we show that 

the mixing between axial and vector GPDs are of order m^/Mjy, which is two orders of 
5^ magnitude less that it follows from the calculation of Mellin moments. 

1 Introduction 

In the last decade many papers were devoted to calculations of the low energy properties of 
generalized parton distributions (GPDs), [BEIEISIEIEIIZIISIEIQIS. The main tool to obtain 
the low-energy expansion of matrix elements is effective field theories based on spontaneous 
symmetry breaking, such as Chiral Perturbation Theory (ChPT) or Heavy Baryon Chiral Per- 
turbation Theory (HBChPT). The pion parton distributions were examined comprehensively, 
but the nucleon GPDs were considered only as Mellin moments. This article is devoted to the 
exploration of the chiral corrections to the nucleon GPDs in the x-space (x is the longitudinal 
hadron momentum fraction carried by a parton). 

Investigation of the chiral structure of nucleon is an important task. Only in this way one 
can obtain the quark mass dependence and the momentum transfer dependence for GPDs 
in a model-independent way. Investigation the chiral expansion gives us the behavior of the 
parton distribution at large impact parameters, where it is governed by pion cloud, see. e.g. 
[TTl rT2] . ChPT provides a good tool for the model- independent investigation of GPD specially 
at x ~ (A x )° (A x ~ M ~ 4:irF n ~ lGeV is the character scale of the chiral expansion). However, 
the most interesting features of the nucleon chiral dynamics lie in the area x ~ m^/A?, ~ 2-10 -2 
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. This region of x is hard-to-reach for ChPT calculations of Mellin moments, but without many 
difficulties can be considered evaluating non-local operators directly. 

GPDs are defined as matrix elements of the light-cone QCD operator. The light-cone separation 
A is dimensional. As a result, one should apply different formal rules of chiral counting for 
different region of x. This leads to ambiguity in restoration of the parton distribution from its 
Mellin moments. Practically, it means that lowering values of x in restored distribution one 
should exponentially increase the number of chiral corrections for the Mellin moments. Working 
directly with non-local operators, one can cover wider areas in x without dramatic increasing of 
difficulties. The price is the incorporation of the non-local operator in the gradient expansion. 

The region of x ~ m^/A 2 ~ 2-10~ 2 is characterized by dominance of the pion cloud contribution 
over the contribution of the hadron core. In the pion case in order to obtain the leading 
chiral correction in this region, one should perform the leading logarithm resummation in 
all orders of the chiral expansion [TUJ [T3|. In the nucleon case, the intermediate scale x ~ 
m 7T /A x ~ 0.15 arises. This scale characterizes the splitting of the pion part of the core from 
the nucleon one. The terms belonging to one or another physical part of distributions have 
well-distinguishable forms in the ^-representation. In the Mellin moment representation, where 
the dimensional parameter A is replaced by the featureless parameter N, it is very difficult to 
distinguish one contribution from another. The Mellin moments project all ranges of x to one 
number. Therefore, extraction of effects at x ~ ^ demands the calculation up to (m n /M) N 
order of chiral expansion. Keeping the light-cone parameter "alive" one does not need such 
extensive calculation, and many results can be obtained already in one-loop. 

The pion-nucleon ChPT has the power-counting problem. The standard tool to avoid the 
counting problem is to use HBChPT, which is a non-relativistic limit of ChPT. Naive attempt 
to set a light-cone operator to the non-relativistic framework is condemned to fail. The infrared 
logarithms of straightforward HBChPT are replaced by some cumbersome functions, which 
often have unnatural singularities in x. The reason is the non-commutativity of the non- 
relativistic limit and the infinite momentum frame nature of the light-cone operator. We have 
consider the non-local operators in the Lorentz invariant pion-nucleon ChPT with subsequent 
reduction to the heavy baryon limit. The heavy baryon limit of the loop- integral is given by 
its infrared part [13]. But loop-integrals with non-local operators are infrared finite even in 
the regime M — > oo, which makes the straightforward calculation wrong. We describe the 
non-contradictory heavy baryon reduction for loop-integral with non-local operators in sec. III. 

In this paper we calculate of the leading chiral correction to nucleon GPDs at £ = directly 
from non-local operators. Also we restrict ourself to leading in the chiral counting, operators. 
The structure of the article is following: in sec. II we derive the effective non-local operators and 
show the difference between our operator and the operators constructed in HBChPT. In sec. 
Ill we perform the detailed investigation of the heavy baryon limit for the non-local operators, 
and derive rules for self-consistent non-local operator consideration within ChPT. The results 
of our calculations are presented in sec. IV. 
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2 Nucleon GPD in the chiral perturbation theory 



Nucleoli GPDs parameterize the matrix elements of non-local operators: 

dK - ixXP+ (m^)rirM-^)\p) = 



2tt 



(1) 



H(x,^A 2 ) + ^^E(x,^A 2 



2M 



t a u{p), 



dX — AP+ (Plg(^)^7 5 A(-^)b) 



2tt 



(2) 



-u{p' 



(nA) 



^H(x,^A 2 )+ l5 ^E(x,^A 2 ) 



T A U( P ), 



where n M is a light-cone vector projecting the "plus" component of momenta; u and u are 
nucleon spinors; M is a nucleon mass. We use the standard notation for kinematical variables 
in the Breit reference frame: 



V 



p + 



A 



p = P 



A 



A_ 



The variable x has the meaning of a momentum fraction carried by a quark with respect to the 
average momentum P of the nucleon. The index A is equal to zero for isospin scalar GPDs, 
and A = 1,2,3 for isospin vector GPDs. 

ChPT is the low-energy effective field theory of QCD. The lowest-order pion-nucleon Lagrangian 
reads: 



£ = $ (ij^ dfi + r M ) - M + yn M 7 M 7 5 ) * + ^ftr [fyE/^f/t + X ] U + x^)] 



(3) 



with M being the nucleon mass and F n ~ 93 MeV being the pion decay constant and x being 
the quark mass matrix. We use the standard notation for the pion field constructions 



u 



U = e F * 



= iu^dfJJu^ . 



Calculations with Lagrangian lead to the power-counting problem, since the nucleon mass 
is of the same order as the normalization scale of ChPT, M ~ A x ~ 1 GeV. The most popular 
solution of the counting problem is provided by HBChPT [15]. The nucleon field is split onto 
soft and hard components: 

where is the nucleon velocity vector, v 2 = 1. One eliminates the component 1-L V using the 
equations of motion and obtains HBChPT Lagrangian. The pion-nucleon part of HBChPT 
Lagrangian reads 

C Nn = Af v (z^(<9 M + r M ) + g a S^) K, (5) 
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where = \ r f(j yv v v is the spin operator, (vS) = 0. 



Right-hand-side of the matrix elements ([HE]) in the heavy baryon limit becomes: 

dK - ixXP+ (m^)nrM-^)\p) = 



2tt 



2 
1 

Pi 



r+X„T A N v ( H + -^e] + ^4An v t a [S + , (SA)]N v 



AM 2 ' 



M 



2 1 H N v S + r A N v + [E + 



H 



1 + 7 



v 2M 2 " 



(6) 



(7) 



where N v is the heavy component of the nucleon spinor: N v (p) = y 1+7 ~2 T = y ^ + it? 5 • 

In the limit M 3> m, the parameter ( ~ ^ < 1. In this article we put £ to zero, by setting 
A + = 0. The effects related to non-zero £ can be of the leading order in M~ l expansion, and 
they will be considered in a separate article. In the limit M — > 00 the GPDs are 

H{x, 0, A 2 ) = q(x, A 2 ), H{x, 0, A 2 ) = Ag(x, A 2 ), 0, A 2 ) = E{x, A 2 ), (8) 



where g and Ag are the parton distribution functions (PDFs) with transverse momentum 
dependance. GPD E can not be obtain in the limit A + = 0. 



2.1 Matching of non-local operators 

In this subsection we construct the effective vector and axial light-cone operators for the nucleon. 
Our procedure is very close to one presented in [21 E] with minor redefinitions. 

The procedure of the matching QCD light-cone operators to the operators in the chiral effective 
theories is well-known. One should introduce the non-local operator with light-cone separation 
A and transformation properties of the initial operator. We start with QCD operators: 

A (\) = g^rfrV-^), (9) 

d A (X) = g(^)^7 5 A(-^). (10) 

These operators transform in mixed representation of the chiral rotations group. It is convenient 
to deal with the left- and right-handed light-cone operators: 

0£(\) = O a (X) + A (A), Of (A) = O a (X) - A (A), 

which transform as the define representation. 
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The fields involved in the chiral Lagrangian transforms under local SUl(2) x SUr(2) as: 

U^RUL\ X ^R X L\ (11) 

^ —> K(u, R, L) Ru K\u,R,L) = K(u,R,L) u L\ (12) 

where R and L are matrices of right and left chiral rotations, and K is a matrix of spinor 
transformations containing field u. 

The lowest order operators contain no derivatives. The most general structures respecting the 
symmetries are: 

o£{\) = 

d(3da \Fl(P,a) *(x)u\x) ft T —u{y)^{y) + F 2 7 (/3, a) V(x)u\x) T —u{y)^{y) 

o£{\) = 

( - T A T A 

d(3da I Fl(p, a) * (x)u(x) i —u\y)^ \y) + F{{p , a) ${x)u{x) $ 7 5 — u\y)^(y) 

where x = Xn^^- and y = \n^^-. The integration area over a and P is \a\ +/3 < 1, these limits 
are provided by additional requirement of polynomiality. The functions F ! (P,a) represent 
generation functions for the tower of low energy constants, and have the meaning of the double 
distributions for the nucleon GPDs with given isospin number in the chiral limit 



The properties of operators fl9lfT0|) under parity transformations demand that F{ = F% and 
F£ = —F[. Thus, the expression for vector and axial light-cone operators are 

A n (X) = J d(3da^(F(((3,a)y(x) jit A (x,y)^(y) + F^{P,a)^(x) ft-ft^(x,y)^(y)) (13) 

OiM = J dpda^(F^P,a)^(x)ji 1 H A (x,y)^(y) + F((P,aMx) ? it A (x,y)^(y))(U) 

where the scalar and pseudo-scalar combinations of pion fields are 

t±(x, y) = u\x)T A u(y) ± u(x)r A u^(y). 

In expressions (I13JI14P one has freedom to choose the normalization of the generating function 
in a convenient way. To recover that the F\ t 2 are the double distributions in the chiral one 
should calculate the matrix element J Q e - txXP + (^p'\0^ n \p) at the tree level, and compare the 
result with definition (JTJ). We obtain: 

dp da F((p, a)5(x - a£ - p) =H (x,^,A 2 ), (15) 

J-l J-1+/3 

where the symbol ° denotes the chiral limit of a quantity. The similar consideration for function 
F 2 (P,a) gives: 

dp da Fi(p,a)5(x-a£- P) =H I (x,^A 2 ). (16) 

1 J-l+B 
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In the limit £ = 0, which is of the special interest in this paper, the a-dependence of the 
operators can be moved out by the translation. The resulting generating functions are 

daF(((3,a) =q I ((3), / daF^(f3,a) =Aq I ((3), (17) 

1+13 J-l+P 
o o 

where Q {(3) and Aq ((3) are the nucleon PDFs in the chiral limit. The normalization of the 
PDFs follows from the vector and axial current operators, it reads: 

d(3°q I ((3) = l, [ df3 Aq ((3) = g a . (18) 
i J -i 

The reduction of the operators (113)141) to the heavy-baryon form is straightforward: 

O^(A) = /d/M^^^ (19) 

+ Fi((3,a) M v {x)S + t A {x,y)M v {y) ) > 



0£ w (\) = I dpdae iM ^ x -^[ V ^Fl{^a)M v {x)t A {x,y)NM (20) 



+ Fl(^a)N v (x)S + t*(x,y)NM\+o(^ 



We remind that in HBChPT nucleons are taken to be non-relativistic, and nucleon momenta 
which flow through the diagrams are reduced momenta: r M = — Mv^. The expressions for 
operators (II 9)20 p have the exponential factor exp(iM(v(y — x))), which restores the reduced 
nucleon momentum to the complete nucleon momentum in the operator vertex. In principal, 
this factor can be absorbed to the definition of the low-energy generating functions F 12 , but this 
would lead to the inconvenient scale of parameters (3 and a. We keep the exponential factors 
as they are, in order to have the simple interpretation ( JT5lrT8|) for the generating functions. 



The matrix element of operators ( ITSirnp contributes to the GPDs H and H at the tree order, 
and to all set of GPDs at one-loop order. The leading contribution to the GPDs E and E is 
given by the operators of the higher order in chiral counting. The operator giving a chiral limit 
of GPD E is proportional to n tJ- d u ('^a tJ-u ^), and it does not give a contribution to H and H at 
tree and one-loop level. Therefore, in this article, which is mostly devoted to GPDs H and H 
we leave such kind of operators out. 

Pure pion non-local operators respecting the transformation properties of ( 19]U0]) can be con- 
structed. Their form and properties were derived at [8] . The pure pion operators for the vector 
and axial PDFs have the form: 

Ot{X) = ^ I d(3da r(P,a)tr [r A (u(x)¥+tf(y) + rf(xff+U(yj)] , (21) 
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OtW = / dfida F(0,a)tr [r A (u{xf$ + U\y) - U\x)H ' + U{y))] , (22) 



^ y y «; 

where the integration area for a and (3 is |a| + < 1, and d = d — d . The axial pion 
operator gives no contribution at A + = (up to one-loop level), and therefore, we will skip it 
in current consideration. 

The low-energy generation function J-(/3, a) belongs to the pion GPD in the chiral limit via the 
expressions: 



o rl rl-P 

H I (x,0= d(3 dadf3J rI (f3,a)[5(x~a^-f3)-(l-I)^5( y x-^( y a + f3))}, (23) 
J-i J-i+p 



where / stands for the isospin, and J 77 " 1 ^/?, a) = ~ (J-(P, a) + (— )J-(— /3, a)). For the forward 

o o 

limit of the pion PDF we will use the notation H 1 (x,£ = 0) =Q I (x). The normalization of 
pion PDFs is: 

dPQ\0) = 8 1 ' 1 . (24) 



2.2 Operators for the Mellin moments 

The Mellin moments of GPDs are widely used in phenomenology The operator of the Mellin 
moment can be found using the well-known relation: 

L ix x " I f/' lx ' Pt « (-t) ^" (t) = ^« (oh+ " d+ Y qioy (25) 

In the usual perturbation theory calculations with a local operator is equivalent to the cal- 
culations with its non-local analogue, they can be compared with each over via direct and 
inverse Mellin transformation. But for the low-energy expansion the situation is not so simple. 
The point is that the left-hand-side of equation (125]) contains the dimensional parameter A, 
which has its own chiral counting. The different areas of A are responsible for the physical 
content of the different regions of x. Dominance of one or another region in A on left-hand-side 
of (1251) leads to different formal counting rules for the derivatives on right-hand-side of (1251) . 
In other words, we can say that the parameter x has it's own counting rule. And applying 
different counting rules for the right-hand-side of (125]) . we recover leading terms responsible 
in corresponding region of x. In this section, we want to carry out an analogy between the 
counting rules and different orders of x, and also shown which part of the nucleon distribution 
was considered in articles [31 SI El El EI - 

The standard counting rules of HBChPT suggests the following: 

<V~0(1), d^~0{q), d^ v ~0{q) (26) 
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where g is the parameter of gradient expansion, i.e. the small momentum or the pion mass 
q ~ m n . We apply these rules to the operator f TTT?]) and obtain in the leading chiral order: 



O 




1 



df3da e 



v+F^a) tf v (0)ti(0)M(0) + 2F 2 (P,a) ST V (0)5+^(0)^(0) 



Taking into account that Mv + = P + and (TTTT) we found 



P+ / ^e-^'O^ m (X) = v + q (x) AT V (0)4(0)K(0) + 2 Aq (x) K(0)S + t A (0)K(0). (27) 



Operator (|27|) allows to obtains only such corrections to PDF those do not touch the intrinsic 
dynamics of the nucleon. In other words using rules (1261) one indirectly uses the approximation 
q(x, A 2 ) = q(x)F(A 2 ), where F(A 2 ) is the form-factor with corresponding quantum numbers. 
This odd result follows from the incorrect counting rules. 

To derive the proper counting rules we start from the statement that in non-local operator 
one can use the standard counting rules for HBChPT. Then using the relation (1251) we can 
obtain the formal counting rules for the Mellin moment operators, which should be used to 
take into account the effects of the original distribution up to certain x. At that, we want to 
keep the standard counting rules of ChPT unchanged, to be able to apply the standard gradient 
expansion. Therefore, we attribute the different counting rules to the light-cone vectors n, h. 
Indeed, after the QCD factorization the light-cone operator does not contain the hard momenta, 
but it still knows about them. In the operators (fl|2|) the dependence on hard momenta is hidden 
in the dependence on the light-cone vector n M . Therefore, the counting rules for "plus", "minus" 
and transverse components of derivative are different for different regions of A. Let us consider 
different regions of x and derive the corresponded counting rules one-by-one. 

The region of x close to unity corresponds to the counting rule n ~ 1. It is equal to statement 
that the short A (A ~ -p- ~ i) dominates in the integral (123|) . The counting rules for the 
Mellin moments to describe this region are the standard counting rules of HBChPT f[26l) . As 
we have shown these rules lead to the factorization of the x-dependance from the operator (1271) . 
It is indeed true for x ~ 1, when one parton is responsible for the whole nucleon. 

The region x ~ -h sets the condition n ~ g -1 , which results in the following counting rules: 



and similar for the nucleon field. Also the consistent counting rules should be attached to the 
vector 



In this region all light-cone derivatives become of the same importance. The operators for Mellin 
moments can be obtained directly by the expansion of the operators (I19ll20p . Numerically this 



<9 + 7r~0(l), <9_tt ~ O (g 2 ) , <9 ± 7r~0(g), 



(28) 




(29) 
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Figure 1: Five diagrams contributing to the leading order. The solid and the dashed lines 
denote proton and pion fields, respectively. 

area corresponds to x ~ ^ ~ 0.15. In this region the loops with composite operators should be 
calculated with great care, because usually ultraviolet logarithms of xP + can become infrared. 
And, indeed, such effects take a place and we discuss them in the next section. 

The region x ~ is governed by the higher derivative terms since n ~ — and <9+7r ~ M . The 
integral (125]) is governed by the huge light-cone distance A ~ ■p-. Therefore, the higher terms 
of the chiral expansion are of the same importance as the leading order and should be taken 
into account simultaneously. The terms belonging to this region have the form of 5-functions 
and their derivatives, see the case of pion [HI [TO] • The resummation of such series in all orders 
of chiral expansion leads to the leading term, which describes the long-range behavior of the 
pion cloud. This mechanism is also known as the "chiral inflation", see |17j . 

We stress, that many problems with intermediate regions can be avoided by calculation directly 
the matrix elements of non-local operators. Keeping the parameter A alive one leaves the 
decision of the interesting x- region until the final result. As a price one has to perform more 
involved calculation. 



3 Non-local operators in the heavy baryon approach 

There are five diagrams contributing to the leading order of PDF chiral expansion, fig.l. The 
calculation of these diagrams is almost straightforward. However, the final result contains some 
uncertainties, related to non-commutativity of the heavy-baryon limit and light-cone definition 
of the operator, which we discuss in this section. 

To show features of the calculation with non-local operators in heavy baryon theories we con- 
sider in details the diagram B. For concreteness we consider vector isovector operator. The 
diagram B is interesting because it mixes the axial and the vector PDFs with each other. Such 
a mixing is believed to be strongly suppressed, because the Mellin moments of this diagram 
gives zero at leading order, due to convolution of the nucleon spin-vector with the nucleon 
velocity. Whereas the calculation with non-local operator gives non-zero result already at the 
leading order. 

The straightforward calculation of loop integral with the Lorentz invariant action ([3]) before 
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the renormalization procedure reads: 
of \ 9aNr a hN 

B{x) = —— r(-l + e) (30) 



x / dp Aq (p) 



i 



where 77 = 1 — 77 and e is a parameter of dimensional regularization: d = 4 — 2e. The integral 
in the brackets gives: 

2M 2 f drjfj [MY + m 2 r]] ~ € S(x - rjfi) = 2(M 2 ) 1 " e | [f + a 2 y] ~ e 9(0 < y < 1), (31) 

where y = § and a = j^. This part is proportional to the large mass, and cannot be considered 
as a correction without additional adaptation. 

The standard approach to the heavy baryon theory posses the rescaling of the nucleon fields 
(J3J) and afterwards expansion in the small ratio a ~ 0.15. In the language of the Feynman 
diagrams it means the partial resummation of the perturbative expansion and it was never 
done explicitly. The bypass way was suggested in [TJ]: the heavy baryon result is provided by 
the infrared part of the loop-integral, whereas the ultraviolet part is to be subtracted into the 
higher order coupling constants. This method works perfectly for usual loop-integrals, but is 
hardly applied to expression ( 1301) . The reason is that the 5-function of the non-local operator 
hardly fixes the low-energy behavior of the loop-integral and does not allow to reach the infrared 
pole at M -)■ 0. 

Let us consider the Mellin moment of the expression f[3"Uj) : 

B N = J dx x N B(x) (32) 
g a b N Nr a hN^. 1 ./(m 2 ) 1 ^ 2M : 



^ + £ > (SB + w I ir) rir ' N [MV + m * n] ") 



where = j dx Aq (x)x N . The infrared part of integral can be found by rescaling of the 
Feynman parameter rj = au and extending the upper limit of integration a^ 1 — > 00: 



(33) 



2M 2 jf 1 d V V v N [MY + m 2 7?] = (m 2 ) 1 ^ (-1 + ^ + 2 )^ + 4 ) Q; 2 + Q ^ 

Therefore, does not contain the infrared logarithm at leading order. 
Applying the same reduction scheme to the integral (13T!) literally we obtain: 

2M 2 I drjfj [M 2 ff + m 2 ri] ~ e S(x - rjp) ->■ 2(M 2 ) 1 ~ e |(y 2 + a 2 y y e 6(y < 1). (34) 
Jo P 

This expression differs from the expression ( 13 ip only by the unnatural Heaviside function, 
which allows the absolute value of variable x to be greater then unity. The Mellin moment 
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for 031]) diverges at infinity, but this divergence is regularized by the dimension regularization. 
Subtracting the pole into the renormalization of coupling constants we comes to expression 
(133]) . At the same time expression does not contain 1/e-pole, and the regularization can 
be removed. But then, the Mellin moment would simply diverge. 

The proper way to reduce the integral (T3TT) to the heavy baryon limit is not to extract the 
infrared part of the loop-integral, but to subtract its ultraviolet part. The results of application 
of these procedures differ in the finite part, which can be fixed by choosing of the proper 
scheme. The ultraviolet part of the loop integral is not changed by the non-locality. Therefore, 
the reduction to the heavy baryon limit for non-local operators consists in the subtraction of the 
pole at M — > oo, with the subsequent renormalization procedure. The result for the expression 
(130]) reads: 

NT a nN f 1 o 

B(x) = — —J dpAq(P) (35) 

This expression is of O (a 2 ) order, in spite of the large mass M in the second term. The Mellin 
moments of it are of order 0{a i ) and coincides with ones obtained from (I33p . 

Note, the expansion of (135]) in a gives asymptotic series because the integral over /3 diverges 
at /3 = x for every term of the expansion. This is a bright example of the non-commutativity 
of the large mass limit with the light-cone formulation of the nucleon structure operators. The 
mixing of vector and axial PDFs, which is believed to be of order a 4 ~ 5 • 1CT 4 , is of order 
a 2 ~ 2 ■ 10~ 2 . In the isoscalar case, where this diagram is enforced by the chiral algebra factor, 
the correction is of order 10% for x ~ 0.8. 



4 Leading chiral correction to nucleon PDF 

We have calculated the leading order corrections to the nucleon parton distributions q ! (x, A 2 ), 
Ag 7 (x,A 2 ) and E ! (x,A 2 ). The results were obtained in the Lorentz invariant theory and 
then reduced to the heavy baryon limit, as it is described in the previous section. We have 
also evaluated Mellin moments of the mentioned above PDFs at next-to-leading chiral order 
independently. Here we present the full set of results depending on the isospin and on the 
Lorentz structure. Since the finite part of the loop-integrals can be absorbed into the higher 
order PDFs we skip it. For the interested reader we present diagram- by-diagram listing in the 
appendix A. 

In order to present the result in more convenient way, we use the following notations: for the 
vector PDFs 

q T (%, A 2 ) = J W + _^_jr^(o<|<i) (36) 
x (imc 1 (j, A 2 ) + aV(/3)AC" + Q\p)Ci A 2 )) , 
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- (sk?£ f"(°<|< 1 ) i^ s ' (f ■ A2 ) + «^ (f • A2 )) • < 37 > 

and for the axial PDF 



Ag 7 (x,A 2 ) = Ag 7 (x) (38) 

1 



+ 



(47rF ff ) 



^/!fK o <i< i )( AV ^'& A2 ) +5 ' ( « Aei (i 



r r - - r 

where 9 (x) and Ag (x) are the nucleon PDFs in the chiral limit, and Q (x) is the pion PDF 
in the chiral limit. The subscript / stands for the isospin. The coefficient functions C 1 , C 1 and 
S 1 are listed below. 

For the vector PDF the functions C 1 are: 

w)/ A 2x 3 9a, f i x 2, m2 o 2 f . ym 2 - A 2 [f] - y - r]y 2 ) 

C(y,A) = — — 5(1 -y)m In— - 3# a / d??— + C(NLog), (39) 

2 V Jo y 2 + a 2 y - ^211(11 - y) 

n i ( A 2n f^, 5 S^\xn \ 2, m2 , 2 r , ym 2 - A 2 (f] -y -r]y 2 ) 

C [y, A) = -[! + — ) 5{l-y)m \n—+g a drj— — (40) 

V 1 J V Jo y 2 + a 2 y - i-FzViV - y) 



M 2 

+£>(NLog). 



For the axial PDF the functions C 1 are 



C°(y,A 2 ) = C%y,A 2 ) + 6g 2 a M 2 y\n(l + ^j+0(NLog), (41) 

C\y,A 2 ) = C 1 (y,A 2 )-2g 2 a M 2 y\n^l + ^fj+0(NLog). (42) 
The parity-mixing functions AC 1 are the same for the axial and vector case: 

AC°(y) = AC (y) = -^6(l-y)m 2 \n%-3g a M 2 y\n(l + ^)+O(NLog), (43) 

2 /i 2 \ y J 

AC\y) = AC\y) = ^6(l-y)m 2 \n% + g a M 2 y\n(l + ^)+0(NLog). (44) 

2 /x 2 V v J 

The pure pion contribution to the vector PDFs is given by: 

Cl{y,A 2 ) = M*(y) f d V R( V )\n^- (45) 
^ Jo Z 1 



/ a 2 v\ P m 2 - nA 2 

+6g 2 a M 2 y In 1 + -J- + / d„ - ^ m2 _^ + O(NLog), 

Cl(y,A 2 ) = (l-g 2 a )5(y) f drjR(rj) In ^ 

Jo H 

-4 9 2MVn (l + J) " *&C*l + £>~'%-*> + ° (NL ° g) - 
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The coefficient functions S denning the vector PDF E(x, A 2 ) are following 

Sl(y,A>) = AglG 1 fdr, *™ ~l {m2 " A ^ + O(NLog), (48) 

where 

T° = 3, T 1 = -1, G° = 3, G 1 = -2, 
i2(»7) = m 2 - 777/A 2 , a = 

The presented results can be compared with existent in literature results for Mellin moments. 
The moments N = of parton distributions q(x, A 2 ) and Aq(x, A 2 ) diagram-by-diagram re- 
produce the vector and axial form factors given in [IB] . The zeroth moment of E(x, A 2 ) differs 
from the corresponding form factor calculated in [18j in logarithm part. Our results reproduce 
the Mellin moments presented in [7] and nucleon gravitational form factors given in pQ, except 
for the contribution of diagram D to E I=0 (x, A 2 ), which differs by factor (—2). Possibly, the 
disagreement in PDF E(x, A 2 ) comes from the non-local operators of higher order, which we 
do not consider in this article. 

The last two terms in (I45H46I) have a special behavior in x. At x ~ 1 these terms contribute to 
the parton distribution q(x, A 2 ) as a 2 , but at x ~ a their contribution increases and becomes 
of order a. In the region x ~ a 2 these terms give contributions of order a 2 In a. The fact 
that corrections change their chiral order with x is a character feature of the pure pion sector. 
Similar terms in nucleon-pion sector (I39ti44|) do not change the order of correction with the 
changing of x. 

Also we are drawing attention to the 5-function contributions in expressions (I45II46I) . These 
contributions have the special nature, because at x ~ a 2 they are formally of the same order as 
the tree term. Together with similar contributions from the higher orders of chiral expansion 
they form the leading chiral correction to the nucleon PDF in the region x < a 2 . The operator 
evidence of the these terms was presented in sec. 2. 2. The resummation of the 5-function 
contributions in all orders of ChPT is an important task. It leads to the model-independent 
picture of the nucleon at large impact parameters. 



5 Summary 

Using nucleon-pion ChPT we have calculated the leading non-analytic chiral corrections for 
nucleon PDFs with non-zero momentum transfer in x-space. We have restricted ourselves 
to the case A + = and to the GPD operators of the leading order. We have derived the 
general rules to deal with non-local operators within the heavy-baryon approach. Our method 
can be applied for operators with various quantum numbers. We have shown that the Mellin 
moments calculated within standard HBChPT approach can be used for restoration of the 
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original parton distribution only in the region x ~ 1, and fails to restore the distribution 
already at x ~ || ~ 0.15. Our results can be used for investigation of the distributions down 
to x ~ ^2 —2-10 2 . We found that the derived corrections give the visible contribution of 
order 10 _1 . 

We have shown that there exist a parity-mixing term, which mixes vector and axial PDFs. The 
mixing contribution is of order a 2 , which can not be deduced from consideration of the Mellin 
moments. Because of uncertainties in the inverse Mellin transformation the mixing terms were 
believed to be of order a 4 , which two orders of magnitude less. 

The pion contributions to the vector GPD contain the term proportional to S(x), which is a 
bright sign of necessity to sum up such terms in all orders of the chiral expansion. They must 
be resummed at x ~ a 2 , since all 5-terms are of the same order as the leading term. The 
resummation can be done in the leading logarithmical order with the help of the recursive 
equations [13] . The result gives the leading asymptotic term in the long distance behavior of 
nucleon GPDs, see e.g.[T7]. 

The slope of PDFs is related to the average impact parameter squared as 

(49) 

A 2 =0 

We have found, that the contribution of the pion sector of PDFs given by (l45H46lT4"8"j) is larger 
in magnitude then the contribution of the nucleon sector (I39H42IT4"T|) for x < a, which is in a 
good agreement with general chiral picture of the nucleon [11]. For example, average squared 
impact parameter (b 2 ^) for the vector isovector parton distribution q I=1 (x, A 2 ) is: 

(b\) = j dx(b\(x)) = -0.056fm 2 + 0.238fm 2 = 0.182fm 2 , (50) 

where the first term comes from nucleon C x {x, A) and the second one from the pion Cl(x, A). 

The important in phenomenology case of non-zero A + , as well as complete calculation of GPDs 
E and E, also can be considered in the presented framework and will be covered in further 
publication. 

We would like to thank M.Polyakov, N.Kivel and G.Gegelia for numerous stimulating discus- 
sions and real interest to our work. This work is supported by BMBF grant 06BO9012. 



A List of diagram contributions 

In this appendix we present the listing of the diagram contributions, in ^-representation and 
their Mellin moments at the leading logarithmical order. The diagrams are shown in fig.l. To 
obtain the results (I36H46I) the renormalization of the nucleon field should be added. 



(b 2 ± (x)) 



dq(x, A^ 
dA 2 
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Diagrams A, B, C, D, E corresponding to the matrix element (jTJ) can be presented as the fol- 
lowing convolution: 



where Q T H (x) and Q E (x) belong to the nucleon parton distributions g 7 (x, A 2 ) and E*(x, A 2 ), 

o o 

respectively. The vector PDF Q (p) should be replaced by the axial PDF Aq I {(3) for the 

o 

diagram B, and by the pion PDF Q\[i) for the diagrams D and E. 
Functions for different diagrams read: 

eUv) = & { *<i-»)=£b4- /%,, ^-f(''; y : ^A +o(NL0g) , 

\ 4 /i 2 Jo y2 + a 2 y _ _ y ) J 

&hAv) = -^(^-^)^l^ + M 2 yln(l + ^)+0(NLog), 

9 

777 

©fr,c(y) = -^(l-yKln — + 0(NLog), 

/i 



</ -+- M2 y M2 



0^(2/) = d^iy) [ dr)R(r])\n^+0(NLog), 

Jo Z 1 

-R(//) = m 2 — r]f]A 2 , a 



where 



M' 

T° = 3, T 1 = -1, G° = 3, G 1 = -2. 

Only triangle diagrams A and D contribute to the distribution E(x, A 2 

rv o,™ 2 _ A 2 



&bAv) = ^ [ V dr ] _ 2 ym 2 A £ y) + 0(NLog) 
Jo r + « 2 2/ - ij^W - y) 

©UlO = /% /^ ) ^ (m2 - A ^ ) + (9(NLog). 



m 2 y m 2 



Three diagrams A, B, C contributing to the axial matrix element (j2J) will be presented in a 
similar convolution: 

Diagiam . 5M|^M £ *, (o £ | < i) AV(«ei, (|) , (52 ) 
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where 6#(V) belongs to the nucleon parton distribution Aq J (x, A 2 ). The axial PDF Ag 7 (/3) 
should be replaced by the vector PDF <? 7 (/3) for the diagram B. 

§i» = ^'(^(l- S )^l 11 ^-m 2 ? j^3 ; + 2M^ln(l + ^)) + 0(NLog), 
eLM = -^(^(l-!/)^lB^ + M 2 S ln(l + ^) + O(NL0g), 



9 

- — ■ 171 

8i, iC (y) = -^(l-y)m 2 ln— + O(NLog) 



Below we present non-analytic contributions to Mellin moments of q(x,A 2 ), Aq(x,A 2 ) and 
E(x, A 2 ) nucleon parton distributions at order 0(a 2 ) in the chiral expansion. We define Mellin 
transformation of the functions Q I H (y), O^iy) and 0#(y) as 

Q I H {N)= f dyy N Q I H {y), Q^N) = f dyy N Q^y), & H (N) = f dyy^y), (53) 
Jo Jo Jo 

and give a list of their moments diagram-by-diagram: 

4 ji z 

9 

ID 

®h,c( N ) = -5 1 > 1 m 2 \n — + 0(a*), 

fx 

l 

&hA N ) = ^Vo^ 2 /^((m 2 -AVl + ^))ln^ + 2m 2 ln^+0(a 4 ), 

o 

®h,e( N ) = ^n,o T^^ln^M + ^a 4 ), 

Jo A 1 

0^(iV) = ^TWln^ + 0(« 4 ), 
0^(iV) = 0(a 4 ), 

2 

- — ■ 771 

e T HC (N) = -c^'Wln— + C(a 4 ). 

/X 

©i;,AW = -^Wln- + 0(« 4 ), 

i 

Q'e^N) = -4^% i0 |^^Mv^+(m 2 -A^)ln^^ 

o 

i 

/ ^ J R(7 7 )ln^M + 0(a 4 ), 
o 
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Functions 0#(iV = 0) and Qjj(N = 0) for zeroth Mellin moment diagram- by-diagram repro- 
duce results of [IS] for the vector and axial form factors. Functions Q^N = 0) differ in the 
logarithm parts. 
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